Coupling of spin and lattice modes in the S=l/2 two-dimensional antiferromagnet 

K 2 V 3 Og with magneto-dielectric couplings 
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Lattice dynamics and magnetic excitations are investigated to elucidate the origin of magneto- 
dielectric effects in the S=l/2 two-dimensional quantum spin compound K2V3O8. We find evidence 
for lattice instabilities at 110 K and 60 K as optical phonon anomalies and a soft mode at 26 cm -1 
in Ai symmetry. Two-magnon excitations in Bi symmetry show an unconventional double-peak 
structure and temperature dependence. This suggests the existence of a split mode near the zone 
boundary caused by a mixing of spin and lattice modes. 



I. INTRODUCTION 

In recent years, the cross-coupling effects of mag- 
netic and electric properties in a single material have 
been intensively investigated due to the fascinating un- 
derlying physics and the potential to realize multifunc- 
tional devices^ Prominent examples are found in multi- 
ferroic materials such as charge ordered LuFe2 04, mag- 
netically driven RMn 2 05 (R=rare earth), and lone pair 
multiferroic compounds R(Fe,Mn)03. 2 ~ 4 In such mate- 
rial classes, low-energy excitations are given by so-called 
electromagnons, which are coupled spin and lattice ex- 
citations. 5 However, it turns out that electromagnons 
are not a monopoly of the multiferroic materials. Elec- 
tromagnons have also been observed in the paraelectric 
phase of a conical-spin magnetically ordered hexaferrite 
Ba 2 Mg 2 Fei 2 22 £ 

Here the question is to what extent phonon-magnetic 
coupled modes retain their hybrid nature in nonmultifcr- 
roic materials but with magnetodielectric couplings. For 
such less strongly coupled materials, magnetism is not 
mutually coupled to ferroelectricity but it is tied to di- 
electric properties. The A 2 B 3 08 fresnoites (A=K,Rb,Cs, 
NH4; B=Ti, V, Mn, Cu) are suitable for studying the 
coupling of a soft phonon mode to magnetic excitations 
because they are susceptible to displace structural phase 
transition and show novel magnetic properties and non- 
linear optical properties^— Actually, magnetodielectric 
effects have been reported in K 2 V30s and are attributed 
to the interplay between structural and electronic degrees 
of freedom^ 

K 2 V30g is the S=l/2 square lattice antiferromagnet 
which has a space group P^bm and lattice parameters 
a=8.870 A and c=5.215 A at room temperatureji 3 - The 
vanadium layer consists of slabs of corner-sharing VO5 
square pyramids and VO4 tetrahedra, separated by non- 
magnetic, interstitial K + ionsJ^ The magnetic behav- 
ior is described by a Heisenberg model with an ex- 
change coupling constant of J = 12.8 K together with 
a Dzyaloshinsky-Moriya (DM) interaction and an addi- 
tional c-axis anisotropy^ The exchange constant is small 



because magnetic (S=l/2) V 4+ -C>5 pyramids are sepa- 
rated by nonmagnetic (S=0) V 5+ -04 tetrahedra. Inter- 
layer couplings and anisotropics lead to an antiferromag- 
netic ordering at Tjy = 4 K, showing unusual field in- 
duced spin reorientations. Magnetization and neutron 
diffraction studies show continuous spin reorientation in 
a basal-plane magnetic fields Theoretical works suggest 
novel magnetic structures, which consist of chiral helices 
with rotation of staggered magnetization and oscillations 
of the total magnetization! 15 i 16 Inelastic neutron scatter- 
ing studies reveal an unexpected mode splitting of a spin 
wave near the (tt/2,tt/2) zone boundary^ 

Thermal conductivity measurements give a hint of a lo- 
cal symmetry breaking as indicated by anomalies around 
110 K. 10 Optical spectroscopy measurements confirm a 
successive structural transition: a distortion of the api- 
cal oxygen of the VO5 square pyramids at about 110 K 
and a weak basal plane distortion around 60 K<ii& Re- 
cent crystallographic measurements unveil a symmetry 
reduction to space groups P^ 2 6c or PJ^nc at 115 KJ£ 

The most remarkable feature is the observation of mag- 
netodielectric effects by magneto-optical experiments: (i) 
a redshift of the 1 eV excitation with a vibration fine 
structure of the frequency 55 cm -1 , and (ii) a narrowing 
of the 2.3 eV excitation with increasing magnetic field^ 
The field-induced changes in the optical excitations were 
attributed to local distortions of the VO5 square pyra- 
mids with external magnetic field. All these seem to 
suggest that a soft lattice is coupled to a spin system 
but clear spectroscopic evidence is missing. Raman spec- 
troscopy can serve as an experimental choice because it 
is extremely sensitive to both lattice anomalies and low- 
lying spin excitations. 

In this paper, we report on the observation of a soft 
mode in Ai symmetry and the two-magnon excitations 
in Bi symmetry which show a similar energy scale and 
temperature dependence. The two-magnon continuum 
exhibits an unexpected double-peak lineshape. These 
feature are interpreted as resulting from the coupling of 
spin and lattice modes at finite momentum. 
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FIG. 1: (Color online) Raman spectra of K2V3O8 in four 
different (a a ), (cc), (a b ), and (ab) polarizations at 5 K 
and 295 K. 



II. EXPERIMENTAL DETAILS 

Single crystals of K2V3O8 were grown in a platinum 
crucible by the flux method. Samples with typical di- 
mensions of 4 x 4 x 0.1 mm 3 were used for Raman scat- 
tering experiments. The Raman spectra were measured 
in a continuous helium flow optical cryostat by varying 
temperature from 1.5 K to room temperature. All spec- 
tra were taken in a quasi-backscattering geometry with 
the excitation line A = 514.5 of an Ar + laser, a DILOR- 
XY spectrometer, and a nitrogen cooled charge-coupled 
device detector. The laser power of 6 mW was focused 
to a 100 /im diameter spot on the crystal. 



III. RESULTS AND DISCUSSION 

A. Lattice instabilities 

Before discussing Raman spectra we present a factor 
group analysis of the expected phonon modes through 
the structural transition. At room temperature K2V3O8 
has a space group PJ^bm and two formula units (Z=2). 



The factor group analysis yields the following irreducible 
representations: 

r tota; = l2A 1 (a'a',cc,E\\c) + 8A 2 + 6B 1 (a'b') 

+ 105 2 (a6) + 21£;(ac,ki,E||a,E||b). (1) 

Here A2+ 2E modes are acoustic phonon modes. 

For temperatures below 115 K a structural transition 
to space groups P4 2 be or PJ^nc takes place. FortheP-^&c 
space group we have the following T-point vibrational 
modes: 

r totQ , = l5Ai(a'a, cc,E||c) + 15A 2 + 15Bi(aV) 

+ lBB 2 (ab) +36^(oc,6c,E||a,E||b). (2) 

For space group PJ^nc the total irreducible representa- 
tions are given as follows: 

Ttotai = 15Ai(aY,cc,E||c) + 15A 2 + 13fli(aV) 

+ l3B 2 {ab) + 32£(ac,&c,E||a,E||b). (3) 

In order to discriminate the phonon anomalies induced 
by the structural phase transition, the Raman spectra 
were compared at 5 K and 295 K in (a a ), (cc), (a b ), 
and (ab) polarizations, as shown in Fig. 1. Here a = a+b 
and b = a—b are rotated by 45° from the respective a and 
b axis. The spectra for the (a'a') and (cc) polarizations 
consist of Ai symmetry modes. The crossed (a'b') and 
(ab) polarizations correspond to the Bi and B 2 mode, 
respectively. The E modes which contribute interplane, 
crossed polarizations could not be studied due to a very 
thin thickness of the crystal. 

At room temperature we observe a total of the 19Ai + 
8P>i +8B 2 modes. The extra 7 Ai +2 Bi modes might be 
due to either local lattice distortions, which activate the 
silent 7 A 2 modes or a leakage of a selection rule. In con- 
trast, the missing B 2 modes are due to a lack of phonon 
intensity as well as to their overlap with larger intensity 
excitations. At T=3 K we identify 36^i + UB ± + 13B 2 . 
The number of the A\ modes is much larger than the ex- 
pected one for both P4 2 bc and PJ^nc. The observed B\ 
and B 2 modes are not compatible with the space group 
P4nc. This indicates that the real crystal symmetry at 
3 K is lower than P4 2 bc and P4nc due to an additional 
base plane distortion at 60 K. Detailed crystallographic 
characterizations are needed to determine the exact crys- 
tal symmetry for temperatures below 60 K. 

The vibrational modes are assigned in comparison with 
the experimental and theoretical Raman works of the 
vanadium oxides AV2O5 (A=Li,Na,Cs,Mg,Ca), which 
has a similar layered structure consisting of VO5 square 
pyramids or VO4 tetrahedra^2r— The high-frequency 
modes at u> > 500 cm -1 represent the V-0 stretching vi- 
brations of the respective VO4 tetrahedra and VO5 pyra- 
mid. The studied compound has two different V sites and 
four different O sites: 01(bridging pyramid and tetrahe- 
dron), 02 (apical to tetrahedron), 03 (bridging tetrahe- 
dra), and 04(apical to pyramid). Since the bond lengths 
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FIG. 2: (a) Temperature dependence of the 891 cm -1 mode. 
Inset: A fit of the 891 cm -1 mode to a sum of two Lorentzian 
curves at 5 K. (b) Temperature dependence of the peak po- 
sition of the 891 cm -1 mode together with a new mode at 
a shoulder, (c) Temperature dependence of the normalized 
integrated intensity of the 891 cm" 1 mode. 



of V-04 and V-02 are 1.582 and 1.628 A, respectively, 
the highest frequency mode of 972 cm -1 observed in the 
(cc) polarization is the V-04 stretching mode of the VO5 
square pyramid. The 926 cm -1 mode is the V-02 stretch- 
ing mode of the VO4 tetrahedra. The 891 and 526 cm -1 
modes are the stretching modes of vanadium and basal- 
plane oxygens. The phonon modes in the frequency range 
between 200 and 500 cm -1 correspond to diverse O-V- 
O and V-O-V bending vibrations. The low-frequency 
u) < 200 cm" 1 are assigned to K atom vibrations and 
torsional and twisting motions of the square pyramids 
and tetrahedra. 

In order to gain local information on lattice distor- 
tions through the structural modulation, we highlight the 
phonon anomalies in the respective low-, intermediate-, 
and high-frequency regimes. Figure 2(a) shows the tem- 
perature dependence of the 891 cm" 1 mode. At room 
temperature the peak consists of a single Lorentzian line. 
Upon cooling through 110 K a new peak shows up at 
the shoulder of the main peak. The temperature de- 
pendence of the frequency exhibits a slight change in 
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FIG. 3: (a) Evolution of the 375 cm -1 peak which appears 
for temperatures below 60 K. (b) Temperature dependence 
of the normalized integrated intensity of the 375 cm -1 peak. 
The dashed line is guide for the eyes. 



slope around 110 K and a weak jump in the vicinity 
of 60 K. The normalized integrated intensity also shows 
anomalies at the respective temperature. The phonon 
frequency anomalies observed below 110 K are correlated 
with the dramatic contraction of the a-lattice parameter 
at the respective temperatures^ Here we note that the 
infrared spectroscopy also observed the splitting of the V- 
04 stretching mode for temperatures below 110 K, which 
was ascribed to the displacement of the apical oxygen 
from the apex of the VO5 square pyramid^ 

The Raman scattering intensity is given by the change 
of a crystal polarizability with respect to the displace- 
ments of a normal mode. Dipole matrix elements and 
interband transition energies are two decisive factors in 
determining the intensity.-- Noticeably, the temperature 
dependence of the 891-cm" 1 phonon intensity resembles 
with that of the peak position of the V 4+ <i — > d transition 
[Compare Fig. 6 of Ref. 18 with Fig. 2(c)]. This suggests 
that local distortions of the VO5 square pyramids are 
closely related to electronic properties. 

Figure 3 shows the intermediate-frequency mode at 376 
cm" 1 , which is induced by the structural modulations for 
temperatures below 60 K. Since the new phonon mode 
corresponds to O-V-0 bending vibrations, it results from 
the basal plane distortions. With decreasing temperature 
the intensity increases almost linearly without saturation 
[see Fig. 3(b)]. This indicates that the base plane mod- 
ulations grow continuously all the way down to the base 
temperature. 

The most striking lattice anomalies are observed in 
the low-energy phonon modes in (cc) polarization as dis- 
played in Fig. 4. The 95 cm -1 mode shows a splitting 
with decreasing temperature through 110 K. It is remark- 
able that the infrared-active mode at 55 cm -1 also ex- 
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FIG. 4: (Color online) (a) Temperature dependence of the 
low-frequency Raman spectra in (cc) polarization, (b) Peak 
splitting of the 95 cm -1 mode as a function of temperature. 



hibits a similar splitting at the respective temperature. 12 
Since these low-frequency modes are associated with the 
torsional and twisting motions of the square pyramidal 
and tetrahedral building blocks, the structural modula- 
tions are a mixture of out-of-plane bending and rotations 
of the vanadium oxide polyhedra. We can summarize the 
evolution of the structural modulations as follows. The 
rotational modulations of the vanadium oxide polyhedra 
start around 110 K and then the base plane distortions 
follow at 60 K. The hardening by 4-5 cm -1 indicates that 
the twisting degree of the V pyramids and the tetrahedra 
increases substantially upon cooling from 110 K. 



B. A soft mode and magnetic excitations 

We turn to more direct evidence for a soft lattice. 
In second-order phase transitions there will always be 
a zone-center Ai component of the soft mode below the 
critical temperature^ Indeed, we observe the sharp peak 
at 26 cm -1 in the Ai symmetry as shown in Fig. 5. Upon 
approaching the phase transition temperature of 60 K the 
26 cm -1 mode softens and then evolves into quasielastic 
scattering. We find no evidence for the presence of a 
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FIG. 5: (Color online) Evolution of the 26-cm - soft mode 
in Ai symmetry which is sensitive to the structural transition 
at 60 K. Inset: Temperature dependence of the soft mode 
frequency. The solid line is a fit to a power-law. See the text 
for details. 



Raman-active soft mode related to the structural phase 
transition at 110 K. The observed soft mode is assigned 
to the torsional mode that flexes the VO5 square pyra- 
mids out of the plane, which could not be detected by 
optical spectroscopy.— Here we note that the layer com- 
pound SrCu2(B03) has a soft mode at 62 cm -1 as an 
in-phase motion of almost all ions along an out-of-plane 
direction.— This soft mode is a consequence of a contin- 
uous corrugation of the CU-BO3 layer. This analogy can 
be taken as further evidence for the basal plane distor- 
tions. 

The inset of Fig. 5 displays the temperature depen- 
dence of the peak position of the soft mode in a tem- 
perature range of 1.5 - 45 K. It is well described by a 
power law, i.e. w = A \T C - Tf, with A = 4.7 cm^/K 
, T c = 55 K and (3 = 0.448. The observed critical ex- 
ponent is close to a mean-field value of (3 = 1/2. For 
SrCu2(B03) the soft mode undergoes an additional soft- 
ening toward low temperatures due to spin-phonon cou- 
pling. 20 This anomaly is not self-evident for the case of 
K2V3O8 although the small discrepancy between the ex- 
perimental data and the fitted curve might have an ori- 
gin of spin-phonon coupling. This is not surprising if the 
long V 4+ -0-0-V 4+ superexchange paths are taken into 
account. 

We now switch to the low-energy spin excitation ob- 
served in (a b ) polarization. Figure 6(a) compares the 
magnetic continuum in the Bi symmetry with the soft 
mode in the Ai symmetry. The continuum exhibits a 
double peak structure extending from 10 cm -1 up to 44 
cm -1 . The primary peak (PI) is centered at 26 cm -1 
while the secondary peak (P2) lies at 16 cm -1 . The pri- 
mary peak position coincides with the soft phonon mode. 
The magnetic continuum is ascribed to two-magnon ex- 
citations by considering the energy scale of the spectral 
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FIG. 6: (Color online) (a) Comparison of the soft mode seen 
in the Ai symmetry and the two-magnon continuum observed 
in the Bi symmetry. Inset: Temperature dependence of the 
integrated intensity of two-magnon scattering in a tempera- 
ture scale of J. (b) Temperature dependence of two-magnon 
scattering in the Bi symmetry. 



weight and the selection rule (see below). 

In two-dimensional (2D) antifcrromagnets, two- 
magnon scattering can be probed as a double spin-flip 
process via the exchange light scattering.— A dominant 
contribution comes from a symmetry- allowed Bi geom- 
etry^ This selection rule strictly holds for our case. 
The density of states of a two-magnon spectrum is given 
by p2(u) = J2k — 2wfc), where uik is a one-magnon 
dispersion. From the upper cut-off frequency of the 
two-magnon spectrum, we can determine the magnon 
energy at the zone boundary. The obtained value of 
2.8 meV(« 22.4 cm -1 ) agrees perfectly with the neutron 
scattering result (compare to Fig. 12 of RefJ^). 

In addition, we can estimate the exchange coupling 
constant between vanadium spins by resorting to the up- 
per cut-off relation, that is, 2Z c JzS — 4.632 J, where J 
is the exchange integral, z is the number of nearest neigh- 
boring spins, and Z c = 1.158 is a quantum correction to 
a classical value d 8 ' 29 We obtain J 13. 6K. This is in 
reasonable agreement with the value of J = 12.8 K ex- 
tracted by inelastic neutron scattering. 17 The small over- 



estimation is due to the presence of interlayer interac- 
tions and DM interactions. Here we note that a conven- 
tional way of determining J from the two-magnon con- 
tinuum is based on the peak energy, which corresponds 
to uip ~ 2.7 J. 30 However, this method is not unambigu- 
ously applicable to K2V3O8 since there is no well-defined 
single peak unlike other antiferromangetsi ' 30 

Figure 6(b) shows the temperature dependence of 
the two-magnon scattering. The spectrum persists to 
T=60 K in the form of short-range magnetic fluctuations. 
This temperature amounts to 15 Tjy(= 5 J). The survival 
of the zone boundary magnon to exceptionally high tem- 
peratures in an energy scale of Tjy reflects the fact that 
in 2D S=l/2 quantum spin systems magnetic excitations 
are scaled by an exchange coupling constant J while the 
classical Neel ordering is caused by interplane interac- 
tions. Although the spectral weight of the two-magnon 
profile matches well with the spin wave theory, a careful 
inspection of the continuum reveals some anomalies. 

First, the relative intensity falls off strongly with an 
upturn about J [see the inset of Fig. 6(a)]. This behav- 
ior is contrasted by other antiferromagnets, which show 
an increase of a total scattering intensity or a moderate 
decrease in a high-temperature limit i 27 ' 30 

Second, the lineshape shows the double-peak struc- 
ture. Since two-magnon scattering is dominated by zone 
boundary magnons, it might be due to the zone bound- 
ary mode splitting. Indeed, inelastic neutron scattering 
identifies the split modes at the zone boundary centered 
around 21 cm -1 with an energy difference of 4 cm -1 ^ 
Since Raman spectroscopy probes magnetic excitations 
as the two-magnon continuum, the energy difference be- 
tween the two peaks is expected to double that of the 
neutron split modes. The observed separation energy of 
10 cm -1 is slightly larger than twice that of the neutron 
split modes. Its possible reason is discussed below. 

In order to explain the observed split modes Lums- 
den et al^ have discussed several possible scenarios, for 
example, magnon-phonon interaction, structural distor- 
tions, lattice disorders, multimagnon scatterings, and or- 
bital degrees of freedom. However, no definite answer 
can be given. Based on the two-magnon lineshape and 
temperature dependence, we propose a coupling of a soft 
mode to magnetic excitations as an origin. The coinci- 
dence of the energy scales and temperature dependence 
makes it difficult to differentiate two modes using exper- 
imental techniques different from Raman spectroscopy. 
Thanks to the Raman selection rules, however, we are 
able to separate them from each other. 

Here we remind the reader that there are soft optical 
and acoustic modes in structural phase transitions. The 
former can be observed by Raman spectroscopy while the 
latter can be detected by Brillouin scattering techniques. 
Although the dispersion of the soft mode cannot be di- 
rectly determined by Raman scattering experiments, we 
can sketch a typical profile by exploiting the fact that 
the soft optical mode has quite flat dispersion and has 
a T-point energy of 26 cm -1 . The resulting dispersion 
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FIG. 7: (Color online) A sketch of coupling of a soft mode to 
a spin wave near the zone boundary. The spin wave data and 
their fit to linear spin-wave theory (solid line) are reproduced 
from Fig. 12 of Ref.— . The dashed line represents a typical 
optical dispersion of the soft mode, which is based on the 
observed F-point energy of 26 cm -1 . 

curve is illustrated schematically as the dashed line in 
Fig. 7. The soft phonon is presented together with the 
spin wave obtained by inelastic neutron scattering exper- 
iments (squares symbols and solid line in Fig. 7). Judging 
from the comparable energy scale of the soft mode at the 
(0,0) point and the magnon at the (w/2,ir/2) point, it 
is apparent that the soft mode crosses the magnon near 
the zone boundary. This provides a compelling explana- 
tion on the split modes. Since a mixing of the spin and 
phonon mode occurs near the zone boundary, a coupled 
spin-phonon character will be manifested at finite mo- 
mentum. Therefore, we expect no anomalous feature for 
the T-point soft mode as Fig. 5 shows. However, the line- 
shape and temperature dependence of the two-magnon 
scattering will be different from conventional antifcrro- 
magnets. As Fig. 6 shows, the spectral form and the 
scattering intensity reflect a hybrid nature. 



Last, we detail the anomalous double-peak feature in 
relation to the neutron split modes. The zone boundary 
magnons have energies of 19 and 23 cm -1 . The cor- 
responding two-magnon density of states has the upper 
cut-off energies at 38 and 46 cm -1 . This energy corre- 
sponds to 4.632 Jj 28 i 29 Due to magnon-magnon interac- 
tions the two-magon peak is shifted down to 2.7 J for 2D 
antiferromagnets^ This relation enables us to estimate 
the peak energies of the two-magnon spectrum as 22 and 
27 cm -1 . Compared to the observed peak energies of 
16 cm -1 (P2) and 26 cm -1 (PI), only the lower-energy 
peak shows an appreciable energy difference by 6 cm -1 . 
The discrepancy suggests that there are additional soft- 
ening channels, specific to the lower-energy branch. This 
might be due to the fact that the lower branch of the 
hybrid modes has a more phonon character. 



IV. CONCLUSION 

To conclude, we have presented Raman scattering 
studies of S=l/2 2D antiferromagnet K^VsOg- Owing to 
magnetodielectric couplings, this system shows marked 
anomalies in phonon and magnetic excitations. Succes- 
sive structural phase transitions at 110 K and 60 K are 
confirmed by pronounced phonon anomalies in frequen- 
cies, intensities, and numbers. Further, for temperatures 
below 60 K we observe a soft mode at 26 cm -1 in Ai 
symmetry and two-magnon excitations in Bi symme- 
try. The two-magnon continuum shows an unexpected 
double-peak structure and a pronounced temperature de- 
pendence. This is interpreted in terms of a mixing of spin 
and lattice modes near the zone boundary. Our result 
suggests that for magnetodielectric materials a mixing of 
phonon-magnetic modes occurs at finite momentum. 
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